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We studied the electronic structure of bulk single crystals and epitaxial films of Fe304. Fe 2p core 
level spectra show clear differences between hard x-ray (HAX-) and soft x-ray (SX-) photoemission 
spectroscopy (PES). The bulk-sensitive spectra exhibit temperature (r)-dependence across the Ver¬ 
wey transition, which is missing in the surface-sensitive spectra. By using an extended impurity 
Anderson full-multiplet model, and in contrast to earlier peak assignment, we show that the two 
distinct Fe-species (A- and B-site) and the charge modulation at B-site are responsible for the newly 
found double peaks in the main peak above Tv and its T-dependent evolution. The Fe 2p HAX- 
PES spectra show a clear magnetic circular dichroism (MCD) in the metallic phase of magnetized 
100-nm-thick films. The model calculations also reproduce the MCD and identify the contributions 
from magnetically distinct A- and B-sites. Valence band HAXPES shows finite density of states at 
Ef for the polaronic half-metal with remnant order above Tv, and a clear gap formation below Tv- 
The results indicate that the Verwey transition is driven by changes in the strongly correlated and 
magnetically active B-site electronic states, consistent with resistivity and optical spectra. 

PACS numbers: 71.10.-w, 71.30.-|-h, 79.60.-i 


In spite of an extensive research of nearly 80 years, [HQ 
the microscopic description of the metal-insulator (MI, 
strictly speaking a half-metal to insulator) transition in 
magnetite (Fe 304 ) remains one of the most important 
open problems in strongly correlated electron systems. 
Magnetite is a classic mixed-valence oxide which under¬ 
goes a first order MI phase transition on cooling below 
Ty^l20 K, at which the electrical resistivity abruptly 
increases by 2 orders of magnitude. Q Simultaneously, 
the inverse spinel crystal structure chanras from cubic 
to monoclinic symmetry. Verwey et aL[^ interpreted 
this MI transition as a charge ordering of the Fe^’*' and 
Fe®’*' states on the B-site. However, the original Ver¬ 
wey model was later refuted by several experiments!^ 0. 
Recent structural studies 0-0 of high resolution neutron 
and x-ray powder diffraction data have indicated a small 
charge disproportionation of only ~0.3e between B-site 
Fe cations with 2-|- and 3-1- formal valency. 

On the theoretical front, recent studies have em¬ 
phasized the role of orbital order, multiferroicity and 
half-metal accompanying the charge-ordering Verwey 
transition. [SEl And while the role of on-site Coulomb 
correlations as well as long-range Coulomb correlations 
have been invoked, their relations with the experimen¬ 
tal electronic structure are still not conclusive. Thus, 
despite intensive investigations, fundamental questions 
about the Verwey transition such as the origin of the 


charge-orbital ordering, the role of disorder and residual 
entropy, 0 121 the precise mechanism of electronic charge 


transport and the accurate microscopic interaction driv¬ 
ing the MI transition are still debated. While enormous 
progresses were achieved in resolving the low-T mono¬ 
clinic structure 0, H, [l^, the superstructure associated 
with the charge ordering was identified only recently as 
a trimeron consisting of t 2 g minority spin electrons delo¬ 
calized over three B-Fe sites: a donor Fe^’*' and two neigh¬ 
boring acceptor Fe®’*' in the formal valency picture. 0 
Further, using time-resolved soft x-ray diffraction and op¬ 
tical reflectivity, it was shown that the trimerons become 
mobile across the MI transition in two steps process. 14 


From the spectroscopic viewpoint, two major funda¬ 
mental questions about the electronic states remain unre¬ 
solved. The first question concerns the peak assignment 
of the Fe 2p core level SX- photoemission spectroscopy 
(PFS)[Ii m and x-ray magnetic circular dichroism 
(XMCD) in Fe ^ 2,3 absorption edge. 17| It is well known 
that Fe 2p SX-PFS and XMCD spectra of Fe 3 04 ex¬ 
hibit typical spectral shapes, which are frequently used 
as fingerprints of magnetite. These spectra have been 
often interpreted in terms of three distinct clusters (A- 
Fe®’*', B-Fe^"*", B-Fe®’*'), with the total spectrum being a 
sum of three distinct spectra with the stoichiometric ra¬ 
tio 1:1:1 for the different Fe sites. However, this model 
has been applied even to the high-T phase, where all B- 
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site Fe atoms are crystallographically equivalent. Recent 
resonant x-ray scattering measurements, however, have 
questioned this ionic model, d, [l^ Thus, all B-site Fe 
atoms should be considered equivalent above Ty. An¬ 
other controversial question concerns the electronic na¬ 
ture of the high-T phase. In a valence band (VB) PES 
study, it was found a gap of 70 meV opening below Ty 
in the occupied density of states, [l^ This matches nicely 
with half the total gap of 140 meV observed below Ty in 
optical spectroscopy. [20| The authors also found a finite 
spectral weight at Ep above Ty and concluded on a half- 
metallic state with remnant order. [SlIIl In contrast, 
Park et al. reported a change in the gap just below and 
above Ty, but the gap was still finite above Ty. 2^ The 
authors accordingly concluded that the Verwey transi¬ 
tion was an insulator to semiconductor phase transition. 
Recent PES studies also suggested that the Verwey tran¬ 
sition is a small polaron insulator-semiconductor phase 
although clear phonon side band fea- 


23 


transition, 

tures have never been observed experimentally. This has 
resulted into two different interpretations of the Verwey 
transition as an “insulator to half-metal” or “insulator to 
semiconductor” transition. 

In this study, we have critically investigated the Fe 
2p core level and VB electronic states near Ep and an¬ 
swer the two fundamental issues described above. By 
employing bulk-sensitive HAXPES, we show that Fe 2p 
core level spectra of Fe 304 exhibit an new bulk char¬ 
acter electronic state which shows T-dependence, but is 
strongly suppressed in surface-sensitive spectra. The new 
feature is magnetically active as seen in MCD-HAXPES 
measured from epitaxial Fe 304 magnetized films. Com¬ 
parison with the calculated MCD shows that it originates 
in the B-site 3d electrons associated with the valence dis¬ 
proportions. We also observed a weak spectral weight 
transfer but over a large energy scale(^2 eV > 10 times 
the gap energy of 90 meV) and a finite density of states at 
Ep in the high-T phase, indicative of strong correlations 
and a metallic state, respectively. 

High quality single crystals of magnetite were synthe¬ 
sized by the floating zone method and the stoichiometry 
is Fe 304 ,oo 2 -[ 23-[23 The 100-nm-thick Fe 304 films were 
fabricated by a pulsed laser deposition technique on a 
MgO(lOO) substrate. HAXPES measurements for 

the single crystal Fe 304 were performed in a vacuum of 
1x10“^° Torr at undulator beam line BL29XU, SPring- 
8 using a Scienta R4000-10kV electron analyzer. 3ll l32| 
The total energy resolution, AE was set to ^170 meV. 
SX-PES was performed at BL17SU, with AE ~200 meV. 
All PES measurements were done under normal emission 
geometry to maximize the depth sensitivity. A single 
crystal of Fe 304 was fractured in situ and sample tem¬ 
perature was controlled to ±2K during measurements. 
The MCD-HAXPES experiments for 100-nm-thick films 
were performed at the undulator beamline BL15XU(33j| 
of SPring -8 using the circularly polarized x-rays from the 
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FIG. 1: (Color online) Measured Fe 2p PES spectra of single 
crystal Fe 304 . (a) Comparison between Fe 2p HAXPES and 
SX-PES. The spectra were corrected by subtracting a Shirley- 
type background, (b) Enlargements of the energy range cor¬ 
responding to Fe 2p2/3 line, (c) T-dependent Fe 2p HAXPES 
spectra across the Verwey transition. Bottom panel; Differ¬ 
ence of 300 K and 150 K spectra to the 110 K spectrum. 


helical undulator. Photon energy was set to 5.95 keV.|27| 
The MCD-HAXPES measurements of Fe 304 in a rema¬ 
nent state were performed at room temperature. Total 
energy resolution was set to 240 meV. 

First, we present experimental Fe 2p HAXPES and 
SX-PES spectra of single crystal Fe 304 at 110 K ob¬ 
tained with photon energies of hv=1.9A keV and 1.0 keV 
in Fig. 1(a) and (b). The kinetic energy of the Fe 2p 
level of SX-PES and HAXPES correspond to a probing 
depth of ~7 and ~80 A, respectively. The core level 
spectra were normalized at the feature “a”. The SX- 


PES spectrum agrees well with reported results[l^ ll 6 | 
and shows a weak shoulder at 708.5 eV, which is conven¬ 
tionally assigned to the Fe^"*" ions. In comparison, the 
spectral weight of the lower binding energy peak labeled 
“/?” in the main line is found to be strongly enhanced in 
the bulk sensitive HAXPES spectrum. Its energy posi¬ 
tion is the same as that of the weak shoulder observed 
in SX-PES. Furthermore, the presence of the broad low 
binding energy satellite labeled “ 7 ” at 707 eV is also 
clearly seen. The strong enhancement of the (3 feature 
and observation of the 7 feature indicate the importance 
of the HAXPES measurements. The (3 and 7 features 
show a weak but observable T-dependence across Ty as 
shown in Fig. 1(c). These new features cannot be in¬ 
terpreted in the usual cluster model with three distinct 
species for the cation (A-Fe^+, B-Fe^+, B-Fe^’*') because 
the observed enhancement of the feature (3 indicates the 
spectral weight contribution of the B-Fe^’*' must be in- 
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FIG. 2: (Color online) Calculated T-dependent Fe 2p spectra 
of Fe 304 , together with the experimental HAXPES spectra. 


creased by a factor of two, if the conventional cluster 
model assignment is true. 

In order to address the origin of the observed behav¬ 
ior, we used an impurity Anderson model (AIM) with 
full multiplets and included a coherent state to describe 
the spectra. This model is well-established and was suc¬ 
cessfully used to study the HAXPES spectra for vari¬ 
ous materials. Details of the model have been de¬ 

scribed in previous works. 3^-37, li^ Let us first con¬ 
sider the high-T phase, where all H-sites are equivalent. 
Thus in the calculation, tetrahedral Fe 04 , and octahedral 
FeOe clusters with trigonal distortion corresponding to 
the A-site and H-site were used, respectively. In both A- 
site Fe 04 (Td) and B-site FeOe{D 3 d) clusters, we used, as 
basis states, six configurations: idJ, 3d^C_, 

3(i^C^, and 3dJ CL . The 3(i®C represents the charge 
transfer between Fe 3d and the coherent state on the top 
of valence band, labeled C. An effective coupling param¬ 
eter V*, for describing the interaction strength between 
the Fe 3d and coherent state is introduced, analogous 
to the Fe 3d-0 2p hybridization V. The parameter val¬ 
ues for both A- and B-site clusters are summarized in 
the supplementary text 27[. The total spectrum was ob¬ 
tained by making a superposition of the spectra from the 
above two clusters in the ratio 1 : 2 , because the composi¬ 
tion ratio of A-site:B-site is 1:2. A relative energy shift 
between A- and B-site of 1.4 eV was required, reflecting 
the known difference of Madelung potential.!^ Note that, 
in the present case, all the B-site Fe cations are treated 
as equivalent for the metal phase. 

The calculated spectrum above Ty is shown in 



FIG. 3: (Color online) (a) Fe 2p HAXPES spectra of 100- 
nm-thick films of Fe 304 for the circularly polarized light, (b) 
Polarization components of A- and B-site are displayed. 


Fig. 2(a) together with the experimental 300 K spectrum. 
To clarify the peak assignment, the A-site and B-site 
components are also shown in Fig 2(a). The agreement is 
remarkable and the experimental features are reproduced 
well by the calculation over the whole energy range. The 
sharp peak “/3” at lower binding energy originates from 
the B-site Fe atoms and the feature a is a mixture of 
A-site and B-site components. The screening effect from 
the coherent state near Ep leads to the formation of the 
low energy weak and broad satellite “ 7 ”. 

Next, we consider the low-T insulating phase spec¬ 
trum. Since the precise XRD studies have suggested 
the three different type of valences, [ 1 , three kinds 


of parameter sets (i.e. Fe^'^'^^'®“*'-like, Fe'^^'^“*'-like and 
Pg 2 . 8 ~ 2 . 94 -i- 7 j]^g^ were used for the B-sites. 271 As we 
show later, the spectral change of the intensity near 
Ep is very small. Thus we needed to use finite val¬ 
ues of V* even for the insulating phase, indicating 
the remnant of the coherent state near Ep in insu¬ 
lating phase. The total calculated spectrum is ob¬ 
tained by a linear combination of the four compo¬ 
nent spectra with a relative weight of Fe(A):Fe(B^ "^'^^-®- 
like):Fe(B'^^-^-like):Fe(B^'®'^^®^-like)=l:^:^:^, respec¬ 
tively [Fig. 2(b)]. Again, the agreement with experiment 
is remarkable and our calculations nicely reproduce the 
observed T-dependence as shown in Fig. 2(c). 

Further support to our interpretation can be obtained 
from independent MCD-HAXPES experiments of 100- 
nm-thick Fe 3 04 films using circularly polarized light, 
which enable us to obtain clear evidence of the B-site 
contribution in the main line. Figure 3(a) shows Fe 2p 
HAXPES spectra of 100-nm-thick film Fe 304 (probing 
depth ^60 A). The Fe 2p spectra were measured for 
the magnetization parallel and antiparallel to the photon 
helicity, where +/— signs refer to antiparallel/parallel 
alignments. Figure 3(a) clearly shows that the HAXPES 
I~ spectra is very similar to the single crystal HAXPES 
spectra in Fig. 1 and the feature (3 and 7 exhibit strong 
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FIG. 4: (Color online) (a) Measured T-dependence of the 
VB spectra of single crystal Fe 304 across the Verwey tran¬ 
sition, using hard x-ray (/ii/=7.94 keV). (b) Solid, dashed, 
and dotted-dashed curces represent the fitted curves, (aii?s + 
ao)f(EB,0), and ao^Es — Eo)6{Eb — Eo), respectively. 


MCD. The observed MCD behavior is well-reproduced 
by calculated MCD signal using the same parameter set 
as for the metallic phase. Since the B-site magnetic mo¬ 
ment is aligned antiparallel to the A-site moment, the 
observed MCD allows a complete characterization of the 
magnetically active states by comparing with the calcu¬ 
lated spectral components. In Fig. 3(b), we show the 
polarization components of A- and B-site spectra. As we 
mentioned before, the main peak of Fe 2p3/2 HAXPES 
spectrum consists of A- and B-site spectra, in which each 
spectrum has a spin polarization due to an exchange in¬ 
teraction between 2p core hole spin and 3d spins. In the 
B-site spectrum, the lower binding energy peak at 708.5 
eV is mainly composed by the parallel spin component of 
(Fe2p, Fe3d)=(4,, -1^) as seen in the top panel of Fig. 3(b). 
For the A-site spectrum also, the lower binding energy 
peak is composed of the parallel spin components (t, T) 
(see second panel of Fig. 3(b)). The resultant /+ spec¬ 
trum of A-site (t, ^) and B-site (t, >^) shows a main peak 
at 710 eV and a shoulder at 708.5 eV (see third panel of 
Fig. 3(b)). On the other hand, the resultant I~ spec¬ 
trum of A-site (4-, ^) and B-site (),, J--) has a double peak 
structure seen in the bottom panel of Fig. 3(b). It should 
be further noted that the HAXPES spectra of thinner 
10 -nm-thick films did not show a clear low binding en¬ 
ergy feature /3 in the previous report,!^ but a very weak 
MCD signal consistent with the present study. This indi¬ 
cates the importance of HAXPES for bulk-sensitive mea¬ 
surements as well thickness dependence for the intrinsic 
electronic structure of Ee 304 and resolves the debate re¬ 
garding the core level assignments. In retrospect, it is 
clear that the controversies discussed in the literatures 
concerning the interpretation of the core level SX-PES, 
x-ray absorption and XMCD stems from these issues. 

In the following, we will focus on the VB energy region 


with HAXPES. The T-dependent VB HAXPES spectra 
near Ep of single crystal Ee 304 are shown in Eig. 4(a). 
The VB consists of the dominantly Ee 3d band with two 
features; the minority spin t 2 g band near Ep (0—1.2 eV 
binding energy) and a broader feature between 1.2—2.5 
eV binding energies. The broad O 2p—Fe 4s bands at 
still higher energies are not shown, but we have mea¬ 
sured it to be consistent with earlier work. [ 2 ^ The 110 
K spectrum shows negligible intensity dX Ep, indicative 
of an energy gap of 90 meV. On increasing temperature 
across Ty, we see a clear spectral weight transfer from 
high to lower binding energies within the minority spin 
Og band, indicating a finite DOS at Ep at T=150 K, 
and a slight increase in DOS at r=300 K. The finite 
DOS at Ep could be simulated using a superposition of 
a simple power law DOS with a constant [271 l42l| . We 
confirmed the slight increase at Ep between T=150 K 
and 300 K. The T-dependence within the minority spin 
t 2 g half-metallic band is also consistent with the changes 
seen in the trimerons using XRD, including the changes 
seen above Ty. The T-dependent spectra also show weak 
spectral changes up to an energy scale of 2.5 eV binding 
energy, which is more than 10 times the gap energy and 
suggest the importance of strong electron correlations. 
The results are consistent with optical spectra which in¬ 
dicate a weak Drude peak in the half-metallic phase due 
to free carriers, albeit indicative of diffusive dynamics 
like a polaronic metal, as well as changes over large en¬ 
ergy scales ~2.5—3.0 eV. Furthermore, asymme¬ 

try due to electro-hole pair shake-up (the Doniach-Sunjic 
line shape) also clearly increased in the O Is PES for 
high-T phase, indicating the metallic electronic state. 271 


In conclusion, HAXPES is used to show that Fe 304 
exhibits new bulk character electronic states, which are 
strongly suppressed within ^^lO-nm-thick surface. In or¬ 
der to explain the observed T-dependence and the MCD 
signal, we need to go beyond the conventional interpre¬ 
tation. By using the extended AIM, we can consistently 
explain these features in the Fe 2p HAXPES spectra. 
The T-dependence originates in the strongly correlated 
and magnetically active B-Ee 3d electronic states. The 
VB spectra show a gap formation at low-T and spectral 
weight transfer from high to low binding energy side with 
finite intensity at Ep in the high-T phase, confirming the 
half-metallic minority spin t 2 g electronic states. 
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